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(N 1 Abstract 



^ ■ The fragment production in multifragmentation of finite nuclei is affected by 

, the critical temperature of nuclear matter. We show that this temperature can 

| be determined on the basis of the statistical multifragmentation model (SMM) by 

analyzing the evolution of fragment distributions with the excitation energy. This 
method can reveal a decrease of the critical temperature that, e.g., is expected 
!> \ for neutron-rich matter. The influence of isospin on fragment distributions is also 

)Q ■ discussed. 

o 

O PACS numbers: 25.70.Pq, 21.65.+f 

Properties of nuclear matter have been under investigation for several decades (see e.g. 
[0]). Besides their general interest for nuclear physics, these studies are very important for 
our understanding astrophysical objects, such as neutron stars. The information about 
nuclear matter in its ground state and at low temperatures is usually obtained as a 
theoretical extrapolation, based on nuclear models designed to describe the structure of 
real nuclei. 

It is instructive to investigate the thermodynamical properties of neutron-rich mat- 
ter under extreme conditions of low densities and high temperatures. This situation is 
expected, for example, at supernova II explosions and during the formation of neutron 
stars. We believe that the liquid-gas type phase transition is manifested, in this case, in 
the forms of instabilities leading to fragment production. The models used for extracting 
nuclear matter properties in the phase transition region should be capable of describing 
the disintegration of homogeneous matter into fragments. It is also important that they 
should allow to be tested in nuclear reactions leading to the total disintegration of real 
nuclei at high excitation energy. The multifragmentation reactions, which started to be 
investigated experimentally nearly 20 years ago (see reviews 0, [§), suit perfectly for this 
purpose. 

Naturally, a thermodynamical model involved in this kind of analysis has to include 
the ingredients necessary for the description of nuclear matter and to provide a good re- 
production of experimental data. So far, the statistical multifragmentation model (SMM) 
0] satisfies these requirements. The SMM has been designed to describe fragmentation 
and multifragmentation of excited finite nuclei 0, H . It includes a liquid-drop approxima- 
tion for individual fragments which corresponds to the liquid-drop description of nuclear 
matter fl|, [|. This is an essential difference of the SMM from other multifragmentation 
models, e.g. pL H, which do not take the nuclear matter properties explicitly into account. 
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Examples of very successful applications of the SMM for the description of different ex- 
perimental data can be found in [3,9-17]. Furthermore, the descriptions of data with 
statistical models confirm that multifragmentation of nuclei, despite of being a very fast 
process, proceeds under a high degree of thermalization. 

Details of the SMM can be found in || , here we concentrate on parts of the model which 
are important for the following discussion. The model describes the fragment formation 
at a low-density freeze-out (p<l/3po, po ~ 0.15fm -3 is the normal nuclear density), 
where the nuclear liquid and gas phases coexist. The SMM phase diagram has already 
been under intensive investigations (see e.g. Ref. Hl8|). The liquid-drop approximation 
suggests that the fragmentation process is accompanied by an increase of the surface of 
nuclear drops. The surface entropy contributes essentially to the statistical partition sum. 
We should point out that the surface free energy depends on the ratio of the temperature 
T to the critical temperature of nuclear matter T c . In the SMM the surface tension cj(T) 
is given by 



a{T) = (7(0) 



'rp2 _ rp2 \ 5 / 4 



T 2 + T 2 



This formula is obtained as a parameterization of the calculations of thermodynamical 
properties of the interface between two phases (liquid and gas) of nearly symmetric nuclear 
matter, which were performed with the Thomas-Fermi and Hartree-Fock methods by 
using the Skyrme forces ||. In addition the scaling properties in the vicinity of the critical 



point (see ||19|| ) were taken into account. At the critical temperature T c for the liquid-gas 
phase transition, the isotherm in the phase diagram has an inflection point. The surface 
tension vanishes at T c and only the gas phase is possible above this temperature. We 
emphasize that our analysis is based on this general effect and that our conclusions will 
remain qualitatively true in the case of other parametrizations satisfying this condition. 
This surface effect provides an effective way to study the influence of T c on fragment 
production in the multi-fragment decay of hot finite nuclei. 

As was established by numerous studies (see e.g. [3-5,9,15,16]) the mass (charge) dis- 
tribution of fragments produced in the disintegration of nuclei evolves with the excitation 
energy. At low temperatures (T<5 MeV), there is a so-called £7-shape distribution cor- 
responding to partitions with few small fragments and one big residual fragment. This 
distribution looks like a result of an evaporative emission. At high temperatures (T>6 
MeV) the big fragments disappear, and there is an exponential-like fall of the mass distri- 
bution with mass number A. In the transition region T « 5 — 6 MeV, however, there is a 
smooth transformation of the first distribution into the second one. The mass distribution 
of intermediate mass fragments (IMFs, fragments with A = 5 — 40) can be approximated 
by a power law A~ T 0, [5[ p0|| . The r parameter decreases with the temperature, goes 
through the minimum at T m 5 — 6 MeV, and then increases again. The small values of 
t indicate that the probability for survival of the biggest fragment decreases drastically 
with the temperature. This behavior may be associated with a phase transition in finite 
systems. It has been shown in many studies (see e.g. []3], [15], [T7| and references therein), 
that there are numerous peculiarities in this region, such as a plateau-like behavior of the 
caloric curve, large fluctuations of the temperature and of the number of the produced 
fragments, scaling laws for fragment yields, and other phenomena expected for critical be- 
havior. Therefore, the temperature characterizing these phenomena is sometimes called 
the critical temperature for finite systems, and r is considered as one of the critical expo- 
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nents. The SMM can describe the critical behavior observed in the experiments [15-17]. 
However, in the present work we use the r parametrization only for the characterization 
of shapes of the fragment mass distributions. In order to avoid any confusion with the 
standard definition of the critical temperature for nuclear matter, we note, in the follow- 
ing, the temperature corresponding to the critical phenomena as a break-up temperature 
for the disintegration of finite nuclei ||. 

The decrease of the surface energy with increasing T/T c (see formula (1)) influences 
the fragment production and, therefore, can be observed in the fragment distributions. 
In this paper we show that this effect can be used for the evaluation of T c by finding 
the minimum r parameter r m j n and the corresponding temperature T m j n . The physics 
behind the phenomenon is quite transparent: If the contribution of the surface energy is 
rapidly decreasing, a nucleus prefers to disintegrate into small fragments already at low 
temperatures. Simultaneously, fluctuations of size of the fragments increase considerably. 
As a result the mass distribution becomes flatter in the transition region, and this leads 
to a decrease of r. 

The SMM calculations were carried out for the Au nucleus (A =197, Z Q =79) at dif- 
ferent excitation energies and at a freeze-out density of one-third of the normal nuclear 
density. This choice is justified by the previous descriptions of the experimental data 
obtained for peripheral collisions [9,12-16]. Below we present results as a function of both 
the temperature and the excitation energy, since they are related quantities 0. 

We have started by using the standard value of the critical temperature implemented 
in the SMM, T c =18 MeV. This value is consistent with many theoretical studies |6|, EH]. 
In Fig. 1 we show typical mass and IMF charge distributions, (Na) ~ A~ T and (Nz) ~ 
Z~ Tz , at an excitation energy E x =7 MeV/nucleon. One can see from this figure that 
the extracted r and r 2 values are very close to each other, since the neutron-to-proton 
ratio of produced IMFs changes very little within their narrow charge range |5|, As 
seen from Fig. 2, the dependences of these parameters versus excitation energy are nearly 
the same. The parameters obtained for primary hot fragments (excited nuclear matter 
drops) and after their secondary de-excitation (measured cold fragments) are also shown 
in this figure. One can see that the difference between the two cases is smallest around 
the minimum r parameter. Therefore, the r min point is weakly affected by secondary 
processes. 

The critical temperature reflects the properties of nuclear matter, however, these prop- 
erties depend on the composition of this matter. For example, T c tends to decrease for 
neutron rich matter J23J. As it was discussed by many authors, see e.g. pi| , the critical 



temperature can be traced back towards neutron rich matter by studies of the disassembly 
of nuclei far from stability. We consider T c as a free parameter in the SMM and analyze 
how r and the break-up temperature can change. In Fig. 3 we show the results for T c =10 
and 30 MeV. In these cases the evolution of r with the excitation energy is similar to 
the one shown in Fig. 2. However, the values of T m j ra and r m j n are essentially different. 
This reflects a considerable change of masses for the dominating fragments. These values 
are plotted in Fig. 4 versus the critical temperature. It is seen that both parameters 
increase with T c and that they tend to saturate at T c — > oo corresponding to the case of 
the temperature-independent surface. This behavior is expected, since in this case only 
the translational and bulk entropies of fragments, but not the surface entropy, influence 
the probability for the fragment formation. 

In the case of neutron-rich matter, the contribution of the symmetry energy increases 
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considerably. It is necessary to take into account the standard dependence of the r param- 
eter on the isospin of the source, while searching for T c . We performed SMM calculations 
of multifragmentation of 124 Sn and 12A La nuclei, which can be used in experiments [p5fl , 
and compare them with the results obtained for Au nuclei. This Sn nucleus is nearly as 
neutron-rich as the gold nucleus, while the La nucleus is neutron-poor. We have used 
the same model parameters as for the Au case, with the standard T c =18 MeV. It is seen 
from Fig. 5 that the neutron-poor source results in slightly lower microcanonical temper- 
atures in the transition region (E x ~ 3 — 5 MeV/nucleon). In Fig. 5 we show also the 
evolution of the r parameter with the excitation energy. The results for Au and Sn are 
very similar and different from those obtained for the La. One can conclude, that IMF 
distributions approximately scale with the size of the sources, and that they depend on 
the neutron-to-proton (N/Z) ratios of the sources. This is because the symmetry energy 
still dominates over the Coulomb interaction energy for these intermediate-size sources. 
One can see from Fig. 5 that the source with the lower N/Z ratio leads to smaller r 
parameters, i.e. to the flatter fragment distribution. We can explain this as an effect of 
the isospin (i.e. the symmetry energy) on fragment formation: A high N/Z ratio of the 
source favors the production of big clusters, since they have a large isospin. Therefore, 
in the transition region, partitions consisting of small IMFs and a big cluster dominate. 
This leads to a very prominent U-shape distribution with large r. When the N/Z ratio is 
low, the probability for a big cluster to survive is small and the system can disintegrate 
into IMFs, which have a favorable isospin in this case. The dominant fragment partitions 
tend to include IMFs of different sizes, and the fragment distribution is characterized by 
small t. Since the difference in the r parameters between the neutron- rich and neutron- 
poor sources is quite large, it can be easily identified. The evolutions of fragment mass 
distributions caused by decreasing T c and changing N/Z ratio can interfere and, therefore, 
the influence of the critical temperature can be separated only after the comparison of 
experimental data with calculations. 

In view of these theoretical findings it is instructive to demonstrate the possibility 
of the application of such an approach for the analysis of experimental data. Presently, 
there are several experimental analyses aimed at the extraction of the critical exponent 
r in reactions with Au nuclei [T7|, |26|. Those methods are not equivalent to the one 
suggested above, however, they are related to the fragment distributions, and the critical 
exponents can be used for an estimation of r m j n and T c . The extracted break-up excitation 
energies E x vary within 3.8 to 4.5 MeV/nucleon, while r is in the range from 2.12 to 2.18. 
As seen from Figs. 2 and 3, for T c > 18, at these excitation energies r is larger than 
T min- We have performed an interpolation of the SMM calculations for the Au sources 
and found that they fit the experimental values of E x and r if T c is in the range between 
18 and 22 MeV. We have also seen from our analysis that r m j n is lower than the extracted 
r by around 15%. The T c estimated in this way is very close to the standard SMM 
parametrization. This conclusion is supported by the analyses of Refs. [15-17] showing 
that the standard SMM reproduces both the experimental critical exponents and other 
characteristics of produced fragments. It is interesting that in Ref. |T7[ a small critical 



exponent r pa 1.88 is reported for multifragmentation of Kr nuclei, which have N/Z 
ratios lower than Au nuclei. This is an indication of the importance of the isospin effects, 
as discussed above. It is worth noting that recently a very close critical temperature 
(T c ~ 16.6 MeV) was extracted from analysis of the break-up ("limiting") temperatures 
in Ref. |[27|| . It is also in agreement with the temperature T c pa 20 MeV obtained in the 
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experiment of Ref. [|28|] . However, it would be important to identify regular changes of T c , 
which needs involving new sources with different isospins. 

In summary, we have pointed out that within the SMM the critical temperature of 
nuclear matter can influence the fragment production in multifragmentation of nuclei 
through the surface energy. We have suggested that this influence can be observed in the 
A~ T parameterization of the fragment yields by finding the minimum r parameter. In the 
experiments the measured values of the parameters are consistent with the standard SMM 
assumption and slightly higher values of the critical temperature, T c « 18 — 22 MeV. The 
SMM predicts that variations of r m i n are especially large in the region of low T c . Therefore, 
there is a possibility to investigate the decrease of T c for nuclear matter under extreme 
conditions, by studying the evolution of r min and T min in the multifragmentation of finite 
nuclei. This could be realized in the case of the neutron-rich nuclei delivered by current 
accelerators with radioactive heavy- ion beams |25|]. The isospin of the source influences 
also the fragment production through the symmetry energy. We have demonstrated how 
the isospin affects the fragment distributions, that should be taken into account in these 
studies. 
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Figure captions 

Fig.l: Average fragment mass and charge yields (Na) and (Nz), after multifragmen- 
tation of Au nuclei at an excitation energy of 7 MeV/nucleon. Solids lines are ~ A~ T and 
~ Z~ T * fits of the IMF yields. 

Fig. 2: Evolution of r 2 (top panel) and r (bottom panel) parameters with the ex- 
citation energy E x of Au sources calculated with the standard SMM parameterization. 
The open circles are for hot primary fragments and the full squares are for observed cold 
fragments. 

Fig. 3: Evolution of r parameters in the SMM for cold fragments with the excitation 
energy E x . The top panel is for the critical temperature T c =10 MeV, the bottom panel 
for T c =30 MeV. 



Fig. 4: The minimum r m j n parameter for cold fragments (top panel) and the corre- 
sponding temperature T m , n (bottom panel) as function of the critical temperature T c of 
nuclear matter. 

Fig. 5: The temperature (top panel) and r parameters for cold fragments (bottom 
panel) versus excitation energy in multifragmentation. The solid, dashed and dotted lines 
are SMM calculations performed for sources with different sizes or isospin (see the figure) 
with the standard T c =18 MeV. 



6 





2 4 6 8 10 12 14 16 18 20 



Z 




5 10 15 20 

E x (MeV/nucleon) 




E (MeV/nucleon) 




2 4 6 8 10 12 

E x (MeV/nucleon) 



